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We investigate the Holstein model by using the dynamical mean-field theory combined with 
the exact diagonalization method. Below a critical temperature Tcr, a coexistence of the pola- 
ronic and the bipolaronic solutions is found for the same value of the electron-phonon coupling 
g in the range QciiT) < g < gc2{T). In the coexistence region, the system shows a first order 
phase transition from the bipolaronic to the polaronic states as T decreases at T = Tp{< Tcr), 
where the double occupancy and the lattice fluctuation together with the anharmonicity of the 
effective ion potential change discontinuously without any symmetry breaking. The obtained 
bipolaronic transition seems to be consistent with the rattling transition in the /3-pyrochlore 
oxide KOS2O6. 
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The electron-phonon interaction in metallic systems 
has been extensively studied for many years because it 
plays important roles in various materials such as A15 
compounds/-* alkali-doped fullerides,^^ magnesium di- 
boride'^^ and manganites.^^ Recently, another class of 
electron-phonon systems in cage structure compounds 
has attracted much interest, where the ion surrounded by 
an oversized cage shows large amplitude local vibrations 
called rattling.^' ^) The superconducting /3-pyrochlore ox- 
ide KOS2O6 with Tc = 9.6K shows a remarkable first or- 
der phase transition without symmetry change at Tp = 
7.5K which is almost independent of the external mag- 
netic field. ^-' Dahm and Ueda^-' have revealed that rat- 
tling phonons with large anharmonicity play crucial roles 
for the anomalous temperature dependence of the NMR 
relaxation rate and the resistivity observed above Tp. 
Therefore, it is important to elucidate the effect of the 
electron-phonon interaction on the rattling transition at 
Tp^^ on the basis of fundamental models such as the 
Hubbard-Holstein model. 

Recently, the Hubbard-Holstein model has been inten- 
sively investigated using the dynamical mean-field theory 
(DMFT)^"' which is exact in inifinite dimensions and is 
expected to be a good approximation in three dimen- 
sions. At zero temperature, a transition from the corre- 
lated metal to the Mott insulator (Mott transition) takes 
place when the on-site Coulomb interaction U is domi- 
nant, while a transition from the polaronic metal to the 
localized bipolaronic state (bipolaronic transition) takes 
place when the electron-phonon interaction g is domi- 
nant. "'^^'^^■' The Mott transition at half-filling is always 
found to be second order. On the other hand, the bipola- 
ronic transition is found to be first order for large values 
of U, while it is second order for small values of U}^' ^^^ 

As for the finite temperature, we have investigated the 
bipolaronic transition in the Hubbard-Holstein model us- 
ing the DMFT and found a coexistence of the polaronic 
and the bipolaronic solutions for the same value of g in 
the range gci{T) < g < 5c2(7')^'*^ (see also Fig. 1 for 
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U = 0). In the coexistence region, the first-order bipola- 
ronic transition takes place at Tp below a critical temper- 
ature Tcr above which the smooth crossover is observed 
instead of the transition. For T > Tp, the system is in the 
bipolaronic state with large lattice fluctuation, while, for 
T < Tp, it is in the polaronic state with relatively small 
lattice fluctuation. The effective potential for oscillating 
ions is double- well type for T > Tp, while it is single- well 
type for T < Tp. When U decreases, the discontinuities in 
physical quantities at Tp decrease together with decreas- 
ing Tcr and the difference \gci — 5c2|- Especially, in the 
Holstein model with U — 0, both of Tcr and \gci — gc2\ 
become very small resulting in numerical difficulties in 
detecting the coexistence region even for T = 0.^^' In 
this paper, we perform the detailed calculations at low 
temperatures down to 1/10* of the bandwidth to obtain 
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Fig. 1. (Color online) The g — T phase diagram of the polaronic 
and bipolaronic states. In the coexistence region g^i < g < gc2, 
both the polaronic and bipolaronic solutions exist. The dotted 
line shows a possible first order transition temperature Tp and 
the cross indicates a critical point at T = Tcr above which the 
crossover between the polaronic and bipolaronic states is ob- 
served. 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



the phase diagrm of the polaronic and the bipolaronic 
states in the Holstein model as shown in Fig. 1. 

The Holstein model is given by the following Hamilto- 
nian 






2^ ekCk<,Ck<T 
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where cj^^ (cJcr) is a- creation operator for a conduction 
electron with wave vector k (site i) and spin ct, and 
flier — Ci^Cia- hi is a creation operator for a phonon at 
site i, where the lattice displacement operator is given 
by Qi = {hi + bi)/y/2ujo. ek, i^o and g are the energy 
for a conduction electron, the frequency of the Einstein 
phonons and the electron-phonon coupling strength, re- 
spectively. 

To solve the model eq.(l), we use the DMFT^") in 
which the model is mapped onto an effective impurity 
Anderson- Holstein model. ^^^ In the case of a semiellip- 
tic DOS for the bare conduction band with the band- 
width W = 1, p{e) = As/l - Ae^/iT, the local Green's 
function G(iw„) satisfies the following self-consistency 
condition, C/o(icj„)~^ = iw„ — f^ ~ (iy/4)^G(ia;„), where 
/x is the chemical potential and Goii^n) is the bare lo- 
cal Green's function for the effective impurity Anderson- 
Holstein model with 5 = in an effective medium which 
will be determined self-consistently. The effective impu- 
rity Anderson-Holstein model with finite g is solved by 
using the exact diagonalization method for a finite-size 
cluster to obtain G(ia;„) at finite temperature T > 0. 
In the present paper, we use 5 site cluster and the cut- 
off of phonon number is set to be 12. We note that the 
numerical results for 6-site are almost the same as those 
for 5-site, and the numerical results for 15 phonons are 
almost the same as those for 12 phonons. We concentrate 
our attention on the particle-hole symmetric case at half- 
filling with (hi) = {J2cr f^io-) ~ 1' ^^"^ ^^ set wq = 0.1. 

The local lattice fluctuation is defined by (Q^) — 
{{Qi — {Qi))"^)- Fig- 2 shows the square root of the nor- 
malized local lattice fluctuation \/{Q^)/{Q^)o as a func- 
tion of g for T = 0.00175, where {Q'^)a = l/2a;o is the 
value for the zero-point oscillation with g — 0. When g 
increases, (Q^) increases gradually for small g, while it 
does steeply at 5 ^ 0.146, and then finally shows a linear 
increase for large g. Remarkably, a coexistence of two so- 
lutions: solutions with small and large lattice fluctuation 
(Q^), is observed for the same value of g in the range 
ffci < .9 < 3c2 as shown in the inset of Fig. 2. For g < gd 
{9 > 9c2), the solution with large (small) (Q^) disappears 
and the solution with small (large) (Q^) exclusively ex- 
ists. In the coexistence region, the system shows a first 
order phase transition which will be discussed later. 

In Figs. 3 (a)-(f), we plot various physical quantities 
as functions of g around 5ci(2) for several temperatures 
T = 0.001, 0.00175 and 0.003. Fig. 3 (a) shows the g 
dependence of the square root of the normalized local 
lattice fiuctuation \/{Q^)/{Q^)o- At low temperatures 




Fig. 2. (Color online) The square root of the normalized local lat- 
tice fluctuation \/{Q^)/{Q^)o as a function of g for T = 0.00175. 
The inset shows the magnification around the coexistence region 
9cl < g < 9c2 where a coexistence of small (Q^) (closed squares) 
and large (Q^) (open squares) solutions is observed. 



T — 0.001 and 0.00175, we observe the coexistence re- 
gion in the range gci{T) < g < (7c2(T), where the range 
for higher temperature T = 0.00175 is smaller than that 
for lower temperature T — 0.001. On the other hand, 
at high temperature T = 0.003, the coexistence region 
disappears and the system shows a smooth crossover be- 
tween small and large lattice fiuctuationg solutions. 

Fig. 3 (b) shows the g dependence of the double- 
occupancy d = {n-ffii). When g increases, d increases 
gradually for small g from the noninteracting value d — 
0.25 (not shown), while it does steeply at g ^ 5ci(2)i and 
then finally shows a linear increase for large g. Similar to 
(Q^), we observe coexistence of a small d solution (po- 
laronic state) and a large d solution (bipolaronic state) 
in the range gci{T) < g < gc2{T) for low temperatures 
T = 0.001 and 0.00175, while a smooth crossover between 
small and large d solutions is observed at high tempera- 
ture T = 0.003. In the bipolaronic state with large d, the 
local charge fiuctuation, {{ui — {ni))'^) = 2d, is enhanced 
together with the local lattice fiuctuation (Q^), while the 
local moment, {sf) — 3(1 — 2d)/4, is suppressed. 

To discuss the spin and charge properties more di- 
rectly, we calculate the local spin susceptibility, x^°'^ = 
Jq {TrSzi{T)s;,i{0))dT with Szi = (nil- ~ '^ii)/2, and the 
local charge susceptibility, Xc°'^ — Jo {TTni{T)ni{Q)) dr , 
and plot them as functions of g in Figs. 3 (c) and (d). 
When g increases, xi°'^ (Xc°'^) decreases (increases) grad- 
ually for small g (not shown), while it does steeply at 
9 ~ 5ci(2): and then finally shows a linear decrease (in- 
crease) for large g. The enhancement (suppression) of 
Xc°'^ (xlf'^) due to the effect of g is consistent with the 
local charge fluctuation (the local moment) mentioned 
above. We note that, x^°^ in the bipolaronic state in- 
creases with decreasing T in proportion to 1/T at low 
temperature as explicitly shown later. 

To examine the correlation between the lattice (or 
the phonon) and the charge fluctuation in more de- 
tail, we calculate the electron-lattice cross correla- 
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Fig. 3. (Color online) The square root of the normalized local lat- 
tice fluctuation \/(Q'^)/(Q^)o (a), the double occupancy d (b), 
the local spin susceptibility x^°'^ (c), the local charge susceptibil- 
ity Xc°'^ (d), the electron-lattice cross correlation function (Tg_; 
(e) and the electron-phonon density correlation function cr^-p (f ) 
as functions of g around the coexistence region for T = 0.001, 
0.00175 and 0.003. 
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Fig. 4. (Color onhne) T dependence of ^/{Q'^)/{Q^}o, d, x^"" , 
Xc°'^, c!e.—l a-iid Te—p for g = 0.1463. The circles (squares) cor- 
respond to the polaronic (bipolaronic) solutions and the closed 
(open) symbols correspond to the stable (metastable) solutions. 
Tp is a possible value of the first order phase transition temper- 
ature. 



tion function, ^^) a^-i — {{b] +bi){hi— {"hi))), and the 
electron-phonon density correlation function/^^ (Je-p — 
{h\hi{hi — {hi))). As shown in Figs. 3 (e) and (f), both 
g dependence of Oe-i and a^-p are very similar to those 
of the local lattice fluctuation and the double occupancy 
(i. e., the local charge fluctuation) shown in Figs. 3 (a) 
and (b) . Then, the strong correlation between the lattice 
(phonon) and the charge fluctuation due to the effects 
of the electron-phonon coupling g is responsible for the 
bipolaronic state with enhanced lattice and charge fluc- 
tuations. 

For various temperatures, we examine the g depen- 
dence of physical quantities to obtain the coexistence 
region g^i{T) < g < g,2{T) (see Fig. 3 for T = 0.001 
and 0.00175). Fig. 1 shows the g — T phase diagram 
where the critical values gci{T) and gc2{T) are plotted. 
When T increases, gc2 (T) decreases monotonically, while 
9ci{T) is almost independent of T, and then gci{T) and 
9c2{T) coincide with each other at a critical tempera- 
ture Tcr ^ 0.0025. For T > Td-, the coexistence region 
disappears and the system shows a smooth crossover be- 
tween the polaronic and bipolaronic states (see Fig. 3 for 
0.003). 

The previous DMFT studies for T = O"'!^) revealed 
that a second order phase transition from the polaronic 
metal to the localized bipolaronic state takes place with 
increasing g a,t g — gc2- At finite temperature below 



Tcr, we find that a first order phase transition between 
the bipolaronic state (T > Tp) and the polaronic state 
(T < Tp) takes place at a transition temperature Tp in 
the coexistence region, similar to the case with the Mott 
transition observed in the Hubbard model. ^''-' Although 
explicit calculations of the free energy to obtain Tp have 
not been done so far, a possible value of Tp is shown in 
Fig. 1. 

In Figs. 4 (a)-(f), we show the temperature depen- 
dence of the physical quantities which are the same as 
those in Figs. 3 (a)-(f) for a fixed value oi g = 0.1463. A 
coexistence of the polaronic and bipolaronic solutions is 
observed for low temperature T < 0.002, while the bipo- 
laronic solution is exclusively observed for high tempera- 
ture T > 0.002. As shown in the g — T phase diagram in 
Fig. 1, it is expected that the system shows a first order 
phase transition between the bipolaronic state (T > Tp) 
and the polaronic state (T < Tp) at a transition temper- 
ature Tp - 0.00175 shown in Figs. 4 (a)-(f). 

In the polaronic state, all of (Q^), d, ^°^ , (Je-i and 
(Te^p gradually increase with increasing T together with 
increasing thermal excitation of phonons, while x^°'^ de- 
creases. On the other hand, in the bipolaronic state, all 
of (<5^), d, x's°'^, cFe-i and ae-p are almost independent of 
T, where the quantum lattice and charge fluctuations are 
fully enhanced due to the strong coupling effects even for 
r = and then the thermal excitation effects are rela- 
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Fig. 5. (Color online) The effective potential for oscillating ions 
VcsiQ) at g = 0.1463 and T = 0.00175 as a function of the nor- 
malized displacement of the ion Q/Qo with Qo = ^/V'^'^O- Solid 
(dashed) line indicates the result of the polaronic (bipolaronic) 
solution corresponding to the open circle (cross) in the inset. 



tively suppressed. We note that the local charge suscep- 
tibility in the bipolaronic state increases with decreasing 
T and shows a Curie law behavior x^c'^ ^ ^/T at low 
temperature where the bipolarons are almost localized. 
This is similar to the case with the localized spins in the 
Mott insulator, where the local spin susceptibility shows 
a Curie law behavior x!°'^ « 1/T. 

Finally, we calculate the effective potential for oscillat- 
ing ions VcsiQ), which is renormalized due to the effect of 
g, by using a variational wave function as previously done 
in refs. [i^-^^)]. Fig. 5 shows VcsiQ) for the polaronic and 
bipolaronic solutions at the same parameters g = 0.1463 
and T — 0.00175, where the first order phase transition 
between the polaronic and bipolaronic states is expected 
to take place as mentioned before. When g increases, the 
harmonic term in Vcs{Q) decreases while the anharmonic 
terms increase (not shown) , resulting in a largely anhar- 
monic effective potential in the strong coupling regime. 
In the coexistence region, the anharmonicity for the bipo- 
laronic state is larger than that for the polaronic state 
as shown in Fig. 5. Then, the first order bipolaronic 
transition is accompanied by the change in the anhar- 
monicity of the effective potential. It is noted that, in 
the Hubbard- Holstein model, Vcs{Q) for the bipolaronic 
state is found to be strongly anharmonic double- well type 
for large C/,^''-' but the anharmonicity decreases with de- 
creasing U together with decreasing the discontinuities 
in physical quantities at Tp, and then Vcs{Q) becomes 
largely anharmonic single-well type for [/ = as shown 
in Fig. 5. 

In summary, we have investigated the half-filled Hol- 
stein model by using the dynamical mean-field theory 
and found that the system shows a first order phase tran- 
sition from the bipolaronic to the polaronic states as T 
decreases at a transition temperature Tp below a critical 
temperature T^r above which the smooth crossover is ob- 
served instead of the transition. At T = Tp, we observe 
the discontinuous changes in various physical quantities 
such as the double occupancy, the local lattice fluctua- 



tion and the anharmonicity of the effective potential for 
oscillating ions, without any symmetry breaking similar 
to the liquid-gas transition. 

In this paper, we set the phonon frequency ujq = O.IW, 
where the value of the critical temperature is obtained 
as Tcr = 0.00251^. For large ujq » W, the Holstein 
model is known to coincide with the attractive Hubbard 
model with the effective on-site Coulomb interaction 
Ues = — 25^/wo. At half-filling, the attractive Hubbard 
model can be transformed into the repulsive Hubbard 
model, ^°'' where the first-order Mott transition takes 
place below the critical temperature Tcr = O.OISW^.^^^ 
Therefore, the bipolaronic transition for luq ^ W is ex- 
pected to take place below the same critical temperature 
Tcr = O.OISW^, although the physics is different from the 
Mott transition. When ujq decreases, the critical temper- 
ature monotonically decreases from Ta — O.OISW^ for 
oJo = 00 to Tcr — >■ for loq — > 0. At the same time, the 
critical values of the bipolaronic transition gci,2 and the 
discontinuities in physical quantities at Tp decrease with 
decreasing wq- We note that, all of Tc,., gci,2 and the dis- 
continuities at Tp increase with increasing U as shown in 
our previous study for the Hubbard-Holstein model. ^^^ 

The first order Mott transition is observed in several 
materials below a critical temperature Tcr, for example, 
the transition metal oxide V2O3 with Tcr ^ 400K^^) 
and the organic conductor k-(ET)2Cu[N(CN)2]C1 with 
Tcr ^ SSK,^^-* where the values of Tcr are consistent with 
the corresponding values of the conduction bandwidth 
of orders of leV for V2O3 and O.leV for organics. In 
these materials, magnetic ordering transition tempera- 
tures are considered to be suppressed due to the frus- 
tration effect and become lower than Tcr, resulting in 
the direct observations of the Mott transition. Then, we 
may also observe the bipolaronic transition when charge 
ordering transition temperatures are suppressed and be- 
come lower than Tcr m strongly coupled electron-phonon 
systems. A promising candidate is the rattling transition 
in KOS2O6. In fact, the observed first order transition 
temperature Tp = 7.5K is consistent with Tcr ~ lOK 
predicted from the calculation with realistic parameters 
of the bandwidth W ~ SeV^**' and the phonon frequency 
ujQ r^ 30meV. 

When the pressure is increased, g/W (U/W) decreases 
with increasing W, and then the first order transition dis- 
appears as clearly understood from the g — T {U — T) 
phase diagram. Such pressure induced disappearance of 
the first order transition is observed in KOs206^^-' as 
well as in ¥203.^^^ In addition, the effects of substitu- 
tion and/or randomness are known to induce the effec- 
tive pressure^^) and is also expected to responsible for 
the disappearance of the first order transition. Actually, 
the rattling transition in KOS2O6 depends on the qual- 
ity of samples and disappears for low-purity samples.^' 
To be more conclusive, we need further investigation on 
the electron-phonon systems with including the effects of 
cage and/or band structures together with the supercon- 
ductivity. 
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